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Suspended-membrane  19-missing-hole  microcavities  in  triangular  lattice  photonic  crystals  are  numerically 
modeled  by  a  three-dimensional  finite-difference  time-domain  method.  The  resonance  frequencies  and  the 
quality  factors  are  calculated  by  interpolation  of  the  discrete  Fourier  transformation  series  of  the  field  with  a 
Pade  polynomial.  The  numerical  results  are  compared  with  the  photoluminescent  spectra  measured  on  the 
cavity  of  a  nearly  identical  dimension.  The  symmetry  properties  of  the  defect  modes  are  analyzed  with  the 
group  theory,  and  resonance  modes  in  the  photonic-crystal  cavities  are  identified  as  irreducible  representations 
of  the  C6„  point  group.  The  far-held  radiations  of  the  identified  modes  in  the  free  space  are  also  calculated  by 
use  of  a  vector  Green’s  function.  It  is  found  that  the  numerical  results  agree  very  well  with  the  experimental 
measurement  in  various  aspects.  ©  2005  Optical  Society  of  America 
OCIS  codes:  230.3990,  140.3460. 


1.  INTRODUCTION 

Modes  in  photonic-crystal  laser  cavities  can  be  engineered 
in  some  detail  through  patterning  of  the  local  dielectric 
environment  on  a  subwavelength  scale.  To  take  advan¬ 
tage  of  this  ability  to  engineer  the  laser  modes,  one  must 
demonstrate  that  experimentally  observed  modes  agree 
with  theoretical  predictions.  Because  the  first  single¬ 
defect  microcavity1  on  a  triangular-lattice  photonic  crys¬ 
tal  has  a  fairly  low  quality  factor  for  resonance  modes  in 
the  photonic  bandgap,  researchers  have  focused  on  other 
cavity  designs,  including  modified  single-defect2,3  and 
larger-defect  cavities.3,4  Recently,  there  have  been  reports 
of  high-Q  modified  single-defect  microcavities  on  graded 
triangular-lattice  photonic  crystals.5  Unlike  a  traditional 
defect  cavity  usually  with  one  hole  or  multiple  holes  re¬ 
moved  from  the  lattice,  those  cavities  are  generally  cre¬ 
ated  by  an  increase  of  the  hole  radius  of  the  defect  and  a 
careful  arrangement  of  the  hole  radii  of  the  adjacent 
holes.  When  the  defect  hole  radius  over  the  photonic  crys¬ 
tal  lattice  is  changed,  the  dielectric  band  mode  is  pulled 
into  the  photonic-crystal  bandgap.  When  the  hole  radii  of 
the  photonic  crystal  cladding,  is  radially  changed  the  pur¬ 
pose  of  minimizing  the  /s-space  overlap  of  the  field  with 
the  radiation  light  cone  is  fulfilled.  In  this  work,  we  focus 
on  resonant  cavities  formed  when  three  concentric  rows  of 
holes  are  removed.  The  resulting  cavity  is  ~2.5  fin i 
across.  There  are  three  reasons  motivating  our  efforts  on 
these  slightly  larger  cavities.  First,  these  cavities  often 
have  a  larger  optical  confinement  factor  than  the  modified 
single-defect  cavities,  since  the  modified  single-defect 
modes  generally  peak  at  the  locations  where  dielectric 
gain  material  has  been  removed.  Because  the  threshold 
modal  gain  is  inversely  proportional  to  the  product  of  the 
quality  factor  and  confinement  factor, 


2t m  1 


in  which  F  is  the  optical  confinement  factor  indicating  the 
percentage  of  the  field  overlap  with  the  gain  medium,  g  is 
the  threshold  material  gain,  n  is  the  effective  index  for 
the  resonance  mode  at  the  wavelength  \0,  and  Q  is  the 
mode  quality  factor,  it  is  useful  to  consider  the  trade  off 
between  quality  factor  and  confinement  factor.  This  trade¬ 
off  is  especially  true,  since  little  is  known  about  the  spa¬ 
tial  dependency  of  the  gain  near  the  etched  photonic- 
crystal  sidewalls.  Second,  the  manufacturing  tolerance 
for  these  slightly  larger  cavities  is  relaxed.  Finally,  elec¬ 
trically  pumped  photonic-crystal  cavities  may  contain  an 
oxide  aperture  whose  minimum  radius  is  likely  to  be  lim¬ 
ited  to  2-3  pm  by  electrical  resistance. 

However,  a  detailed  numerical  analysis  on  such  multi- 
mode  cavities  can  be  a  difficult  problem,  since  there  can 
be  many  closely  spaced  modes  in  these  larger  cavities.  In 
this  paper,  numerical  methods  are  introduced  for  efficient 
characterization  of  multimoded  photonic-crystal  micro¬ 
cavities.  It  is  shown  that  we  obtained  excellent  agreement 
between  our  numerical  predictions  and  the  observed 
modes.  The  paper  is  organized  as  follows.  Section  2  pre¬ 
sents  the  results  of  a  three-dimensional  (3-D)  finite- 
difference  time-domain  (FDTD)  simulation.  As  a  part  of 
quality-factor  calculation,  a  time-domain  discrete  Fourier 
transform  (DFT)  of  the  magnetic  field  component  Hz  at 
low-symmetry  locations  of  the  cavity  is  interpolated  with 
Pade  polynomials.  The  quality  factors  are  then  calculated 
by  the  full  width  at  the  half  magnitude  of  each  resonance 
peak  in  the  Pade  interpolated  frequency  response.  Experi¬ 
mental  data  in  the  form  of  photoluminescence  (PL)  spec¬ 
tra  from  the  fabricated  devices  are  also  given  and  com- 
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pared  with  the  calculation.  In  Section  3,  the  field 
distributions  for  the  resonant  modes  of  interest  are  ob¬ 
tained  by  application  of  a  digital  filter  during  the  FDTD 
simulation.  The  resonant  modes  are  classified  in  the  C6v 
point-symmetry  group  by  projection  of  the  electrical  and 
magnetic  fields  onto  irreducible  representations  of  the 
group.  The  modes  are  invariant  under  the  set  of  symme¬ 
try  operations  of  the  irreducible  representations  of  the 
relevant  point  group.  Section  4  presents  the  free-space 
far-field  radiation  pattern  of  a  doubly  degenerate  reso¬ 
nant  mode  calculated  with  the  vector  Green’s  function  ap¬ 
proach.  The  result  is  compared  with  the  experimental 
measurement,  which  agrees  very  well  with  the  theoretical 
analysis. 


2.  RESONANCE  MODES  AND  QUALITY 
FACTORS 

The  frequency  response  of  the  cavity  is  calculated  by  a 
Fourier  transform  of  the  time-dependent  field  at  the  low- 
symmetry  locations  recorded  from  a  3-D  FDTD 
simulation.6'7  However,  since  the  dielectric  constants  of 
all  of  the  materials  are  assumed  to  be  real  and  constant  in 
the  simulation,  we  are  limited  to  predicting  the  cold  cav¬ 
ity  behavior  of  the  cavities.  Since  only  a  limited  number  of 
modes  exist  in  such  a  device,  the  frequency  resolution  of 
the  temporal  Fourier  transform,  which  is  given  by 

A  f  M 

—  = - ,  (1.2) 

f  N(a/\0)’ 

in  which  a  is  the  lattice  constant  of  the  photonic  crystal, 
A-o  is  the  free  space  wavelength,  N  is  the  total  number  of 
samples  for  the  Fourier  transform,  and  M  is  the  number 
of  time  steps  it  takes  for  the  field  to  travel  a  distance  of  a 
in  the  free  space,  shows  that  the  frequency  resolution  can 
be  made  to  be  much  smaller  than  the  mode  spacing. 

In  contrast,  the  situation  for  the  quality-factor  calcula¬ 
tion  for  the  multimode  cavity  is  more  complicated.  Calcu¬ 
lation  of  the  quality  factor  of  each  resonance  by  the  ratio 
of  the  power  exiting  the  device  and  the  total  energy  in  the 
cavity1  generally  requires  a  significant  number  of  time 
steps,  because  it  requires  the  use  of  a  time-domain  filter 
narrow  enough  in  the  frequency  domain  so  that  only  the 
desired  mode  is  excited.  For  a  given  filter  window,  the  ra¬ 
tio  of  the  filter  bandwidth  BW  and  its  central  frequency  f 
is  given  by 


BW  -n£M 

- = - ,  (1.3) 

f  N(a/\0) 

where  £  is  the  transition  width  factor  of  equivalent  Kaiser 
window.  For  example,  a  Hanning  window  has  a  £  of  5.01. 
The  spacing  between  modes  for  multimode  microcavities 
is  generally  so  small  that  an  implementation  of  the  sharp 
filter  is  impractical.  The  method  is  also  somewhat  sensi¬ 
tive  to  the  filter  efficiency  and  tends  to  be  less  accurate  as 
the  quality  factor  of  the  mode  increases,  since  the  lost 
power  is  very  small  for  a  high-Q  resonant  mode. 

In  this  paper,  because  of  the  difficulty  with  efficient, 
narrow-bandpass  digital  filter,  the  quality  factor  is  in¬ 
stead  calculated  from  the  ratio  of  full  width  at  the  half 


magnitude  of  the  cavity  resonance  in  the  frequency  do¬ 
main,  Aa>,  to  the  center  frequency  oj0.  The  inherent  limi¬ 
tation  of  this  approach  is  its  insufficient  frequency  reso¬ 
lution;  the  smallest  identifiable  spectral  width  usually 
cannot  resolve  a  quality  factor  higher  than  a  few  hundred. 
Distortion  to  the  spectrum  is  also  introduced,  because  the 
numerical  simulation  terminates  before  the  impulse  re¬ 
sponse  is  fully  evolved.  This  has  the  effect  of  viewing  the 
true  time-domain  response  through  a  rectangular  win¬ 
dow,  which  translates  mathematically  into  the  convolu¬ 
tion  of  the  true  spectrum  with  a  sine  function.  The  convo¬ 
lution  widens  the  peaks  in  the  spectrum,  among  other 
effects.  This  distortion  can  be  reduced  with  an  increase  of 
the  time-response  window,  but  at  a  price  of  a  longer  simu¬ 
lation  time.  Pade  interpolation  addresses  the  problem 
through  extrapolation  of  the  electromagnetic  field  in  the 
time  domain  beyond  the  actual  simulation  window.8  The 
DFT  series  P((ok)  is  interpolated  with  a  Pade  function, 
which  is  the  ratio  of  an  order  I  and  an  order  J  polynomi¬ 
als  Qi  and  Dj. 


P(<»k)  = 


Qii^k) 

Dj((ok) 


2 

i=0 

2  PMy 

i= o 


(1.4) 


Assuming  /30  to  be  unity,  coefficients  a,  and  /3,  in  Eq.  (1.4) 
can  be  solved  with  I+J+l  Fourier  transform  samples, 
and  a  continuous  Pade  function  interpolating  the  Fourier 
spectrum  is  derived.  The  accuracy  of  the  frequency  re¬ 
sponse  is  significantly  improved  beyond  the  frequency 
resolution  of  the  discrete  Fourier  transform.  By  combin¬ 
ing  Pade  approximation  with  the  DFT,  one  can  determine 
the  quality  factors  of  all  modes  in  the  cavity  with  just  one 
FDTD  simulation. 

Figure  1  shows  the  TE2-like  (even)  resonance  mode  fre¬ 
quencies  and  their  quality  factors  for  a  triangular-lattice 
photonic-crystal  defect  microcavity,  where  the  2  direction 
is  the  epitaxial  growth  direction.  The  cavity  is  created  in 
a  dielectric  membrane  suspended  in  the  air.  The  mem¬ 
brane  has  a  refractive  index  n  of  3.4  and  a  thickness  d  of 
0.45  normalized  against  lattice  constant  a.  The  defect  is 
formed  when  the  three  concentric  circles  of  air  holes,  a  to¬ 
tal  of  19  holes  are  filled,  with  dielectric.  All  air  holes  have 
a  normalized  radius  of  r/a  =  0.33  in  the  simulation,  except 
those  immediately  adjacent  to  the  cavity  defect,  which 
have  a  slightly  smaller  radius  ofr/a  =  0.3.  For  the  purpose 
of  comparison,  the  PL  spectrum  of  an  actual  photonic- 
crystal  defect  cavity  pumped  below  threshold  is  also 
shown  in  Fig.  1  as  a  curve.  The  devices  were  fabricated  in 
an  epitaxial  layer  structure  consisting  of  a  20-nm-thick 
InGaAs  etch-stop  layer,  a  60-nm  sacrificing  InP  top  layer 
to  protect  the  quantum  wells  during  the  processing,  and  a 
240-nm-thick  InGaAsP  layer  containing  four  compres- 
sively  strained  InGaAsP  quantum  wells  designed  to  emit 
at  1.5  pm  at  room  temperature.  The  cavity  dimensions 
are  8  pm  by  10  pm,  as  shown  in  the  scanning  electron  mi¬ 
crograph  image  in  the  Fig.  1  inset.  In  all  cases,  there  are 
either  six  or  seven  periods  of  a  triangular-lattice  photonic- 
crystal  cladding  around  the  central  defect.  The  fabricated 
devices  had  a  lattice  constant  of  500  nm  and  an  r/a  value 
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Fig.  1.  Calculated  resonance-mode  frequencies  and  their  qual¬ 
ity  factors  for  a  suspended  membrane  19-missing-hole  photonic- 
crystal  defect  cavity  (bars)  with  a  membrane  thickness  of  d/a 
=  0.45  and  hole  radii  of  r/a- 0.33  except  the  inner  holes  around 
the  defect,  which  have  slightly  smaller  r/a  =  0.3.  It  is  compared 
with  the  measured  photoluminescent  spectrum  of  a  cavity  with 
similar  dimensions  (curve).  The  left  and  bottom  axes  are  for  the 
calculation,  whereas  the  right  and  top  axes  are  for  the  measure¬ 
ment  curve.  The  difference  of  frequency  is  within  0.0035  in  nor¬ 
malized  scale,  or  —16  nm  in  free-space  wavelength.  Resonance 
modes  with  a  quality  factor  over  1000  in  the  material  gain  region 
are  also  labeled  with  their  corresponding  irreducible  representa¬ 
tion  in  the  C6„  point  symmetry  group.  The  inset  shows  the  top 
view  of  the  cavity. 


of  approximately  0.33.  The  left  and  bottom  axes  are  for 
the  calculated  quality  factors  and  the  resonance  frequen¬ 
cies,  whereas  the  right  and  top  axes  are  for  the  measured 
optical  intensity  and  frequencies.  The  horizontal  axes  of 
the  two  plots  are  intentionally  shifted  by  less  than  1%,  or 
16  nm  in  wavelength.  The  peaks  in  the  PL  spectrum  gen¬ 
erally  align  very  closely  with  the  predicted  resonances, 
except  for  the  ones  around  the  normalized  frequency 
a/A.0  =  0.33,  which  will  be  explained  in  Section  4.  This  fre¬ 
quency  shift  is  attributed  to  the  fabrication  variation  of 
the  microcavity.  The  case  in  which  two  resonance  modes 
are  measured  for  the  wavelength  in  which  only  one  reso¬ 
nance  is  calculated  is  also  explained  in  Section  4.  It  is 
worth  mentioning  that  the  lasing  mode  of  the  microcavity 
is  the  mode  with  the  highest  predicted  Q  at  the  normal¬ 
ized  frequency  of  a/\0  =  0-304. 

Notice  that  the  PL  measurement  of  the  cavity  is  over¬ 
laid  with  the  calculated  quality  factor  of  the  defect  modes 
instead  of  the  DFT  of  the  time  evolution  of  the  field.  Be¬ 
cause  the  intensity  of  the  Fourier  spectrum  does  not  en¬ 
tirely  reflect  the  photon  lifetime  of  the  resonant  modes, 
the  excitation  efficiency  of  the  initial  condition  with  each 
defect  mode  strongly  affects  the  magnitude  of  the  Fourier 
intensity.  Even  though  high-Q  resonance  has  a  slower  en¬ 
ergy  decay  in  time,  which  results  in  a  higher  Fourier  in¬ 
tensity  at  the  resonance  frequency  given  the  initial  en¬ 
ergy  is  identical,  the  overlap  of  the  random  initial 
condition  might  have  strongly  favored  low-Q  modes.  Ad¬ 
ditionally,  the  amplitude  of  the  peaks  in  the  PL  measure¬ 
ment  reflects  the  out-of-plane  loss,  not  the  total  loss.  The 
predicted  high-Q  modes  are  identified  in  the  optical  spec¬ 
trum  as  sharper  resonances,  not  the  ones  with  higher  op¬ 
tical  intensity. 


3.  CLASSIFICATION  OF  RESONANCE 
MODES  WITH  GROUP  THEORY 

As  can  be  seen  from  Fig.  1,  there  are  usually  tens  of  reso¬ 
nant  modes  existing  for  a  multiple-defect  photonic-crystal 
microcavity.  Not  all  of  the  modes  in  the  photonic  bandgap 
are  interesting.  For  those  that  do  raise  interest,  their  field 
distribution  can  be  obtained  with  a  narrow-bandwidth 
digital  filter.  The  field  distribution  of  a  resonant  mode 
provides  insight  into  its  loss  mechanism.  The  classifica¬ 
tion  of  modes  with  the  group  theory  provides  much 
broader  and  more  general  information,  such  as  degen¬ 
eracy. 

The  TEz-like  (even)  resonant  modes  in  a  19-missing- 
hole  suspended-membrane  defect  cavity  belong  to  the  ir¬ 
reducible  representations  of  the  C6v  point  group.9’10  Table 
1  shows  the  character  table  for  the  C6„  point-symmetry 
group.  Together  with  the  identity  operation  E  that  keeps 
the  field  as  is,  the  C%v  group  is  composed  of  the  following 
symmetry  operations11: 

C6„  =  {E,C2,2C3,2C6,3ad,3av}-  (1.5) 

The  C6,,  point-symmetry  group  consists  of  six  irreducible 
representations  labeled  as  A1;  A2,  Bt,  B2,  Eu  and  E2.  A 
resonant  mode  in  a  19-missing-hole  photonic-crystal  de¬ 
fect  cavity  assumes  the  property  of  one  of  the  irreducible 
representations.  For  TE2-like  (even)  resonance  modes  in 
the  photonic-crystal  microcavity,  the  magnetic  fields  are 
nearly  scalar.  In  particular,  the  magnetic  fields  at  the 
midplane  of  the  membrane  consist  only  of  the  Hz  compo¬ 
nent.  The  vector  electric  field,  however,  is  dominated  by 
the  Ex  and  Ey  components.  The  symmetry  operations  are 
therefore  applied  to  the  scalar  magnetic  field  Hz  and  the 
vector  electrical  field.  Figure  2(a)  shows  the  Hz  compo¬ 
nent  of  the  lasing  mode  at  the  midplane  of  the  membrane 
calculated  by  a  3-D  FDTD  method  with  a  60,000  time- 
step  Blackman12  filter.  Projecting  this  mode  onto  the  six 
irreducible  representations  of  the  C6v  point  group  by  ap¬ 
plying  mirror  and  rotation  operations11'13  to  the  mode 
profile  numerically,  we  obtained  the  projections  of  the  sca¬ 
lar  magnetic  field  component  Hz  onto  the  six  irreducible 
representations  of  the  C6„  point  group.  The  projection  op¬ 
erators  onto  the  irreducible  representations  are  given  by 

pu)  =  ^  x(i\R)Pr,  (1.6) 

"  R 

in  which  PR  is  the  operator  corresponding  to  coordinates 
transformation  R,  lj  is  the  dimension  of  the  jth  irreducible 
representation,  and  h  is  the  total  number  of  elements  in 


Table  1.  Character  Table  for  the  C6v 
Point-Symmetry  Group 
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2C6 
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1 

1 

1 

1 

-1 

-1 

1 

-1 

1 

-1 

-1 

1 

B2 

1 

-1 

1 

-1 

1 

-1 

Ei 

2 

-2 

-1 

1 

0 

0 

E2 

2 

2 

-1 

-1 

0 
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Fig.  2.  Magnitude  of  the  Hz  component  at  the  midplane  for  the 
19-missing-hole  suspended-membrane  photonic-crystal  defect- 
cavity  mode  considered  (top),  with  a  two-dimensional  dielectric 
contour  overlaid.  Its  projections  onto  the  irreducible  representa¬ 
tions  of  the  C6„  group  are  labeled  by  A1;  A2,  Blt  B2 ,  Elt  and  E2. 


nearest-neighbor  approach  for  the  values  that  is  not  di¬ 
rectly  given  by  the  FDTD  simulation. 

The  remaining  projection  also  proves  a  point  we  made 
in  Section  2.  The  quality  factor  of  a  multimoded  cavity 
calculated  with  the  ratio  of  dissipated  power  and  total  en¬ 
ergy  is  distorted  by  the  remaining  energy  of  the  neighbor¬ 
ing  resonances  inadequately  suppressed  by  the  time- 
domain  filter.  Even  with  only  a  fraction  of  a  percent  of 
energy  remaining  in  other  modes,  the  calculation  of  the 
quality  factor  cannot  be  accurate  for  modes  having  a  qual¬ 
ity  factor  larger  than  a  few  thousand.  As  the  filter  band¬ 
width  and  the  sideband  suppression  ratio  generally  de¬ 
crease  inversely  with  the  filter  window  length,  given  a 
filter  window  function,  realization  of  a  desirable  filter  is 
computationally  cumbersome.  Certainly,  one  can  use  the 
result  after  the  projection  operation  as  an  initial  condition 
for  elimination  of  the  noise  from  the  undesired  modes  and 
acceleration  of  the  convergence. 

An  additional  error  occurs  in  the  projection  of  the  vec¬ 
tor  electric  field  as  a  result  of  the  noncollocation  grid 
scheme  employed  in  the  FDTD  simulation.  Because  the 
electric  field  components  Ex  and  Ey  are  sampled  at  differ¬ 
ent  spatial  locations,  interpolation  is  required  for  one  of 
the  field  components  to  obtain  its  value  on  the  locations 
where  the  other  component  is  sampled.  For  simplicity,  a 
linear  interpolation  is  applied  in  the  transformation. 

Accordingly,  the  electric  fields  Ex  and  Ey  of  the  same 
mode  assume  the  symmetry  properties  of  the  A2  irreduc¬ 
ible  representation  of  the  C6v  group,  as  shown  in  Fig.  3. 
Here,  we  followed  the  projection  rules  elaborated  by 
Sakoda,13  that  is,  the  projections  of  the  real  vectors  and 
pseudovectors  conform  to  the  same  formula.  The  differ¬ 
ence  between  the  E  and  H  fields  originates  from  the  fact 
that  the  electric  field  is  a  true  vector,  whereas  the  mag¬ 
netic  field  is  an  axial  vector.  Under  mirror  reflections  or 
improper  rotations,  the  electric  field  behaves  as 

E  — >  E'  =  -  E,  (1.7) 


the  group.  The  projections  of  the  calculated  field  for  the 
mode  that  lases  experimentally  are  plotted  in  the  same 
linear  gray  scale  to  emphasize  the  fact  that  the  fields 
have  a  trivial  amount  of  overlap  with  representations 
other  than  Aj.  Energy  in  projections  other  than  A1  ac¬ 
count  for  less  than  1%  of  the  total.  Additionally,  the  mag¬ 
netic  field  is  nearly  identical  before  and  after  a  set  of  pro¬ 
jection  operations  characterizing  Aj  irreducible 
representations.  Therefore  the  magnetic  field  Hz  of  the 
lasing  mode  is  identified  as  an  A}  mode  of  the  C6v  point 
group. 

The  remaining  projections  are  attributed  to  the  ineffi¬ 
ciency  of  the  temporal  bandpass  filter  and  the  numerical 
errors  involved  with  the  symmetry  operations.  Even 
though  a  Blackman  filter  has  a  sidelobe  suppression  ratio 
of  57  dB,  some  traceable  amount  of  energy  still  remains 
after  the  bandpass  filtering,  depending  on  the  initial  ex¬ 
citation  condition  and  the  photon  lifetime  of  adjacent 
resonant  modes.  In  addition,  the  rotations  and  the  mirror 
reflections  of  the  multiples  of  60  deg  cannot  be  imple¬ 
mented  in  the  Cartesian  coordinates  without  the  applica¬ 
tion  of  antialiasing.  Here,  we  have  applied  a  weighted 


whereas  magnetic  field  behaves  like 

H^H'  =  -H.  (1.8) 

The  characters  for  irreducible  representations  Aj  and  A2 
have  opposite  signs  for  mirror  reflections  crd  and  cr„.  Some 
other  literature  applies  an  additional  negative  sign  to  all 
improper  rotation  operations  of  the  pseudovectors.  The 
electric  and  magnetic  fields  therefore  conform  to  same  set 
of  irreducible  representation.  In  this  paper,  we  will  label 
the  resonant  mode  based  on  the  representation  of  its  mag¬ 
netic  field.  For  resonant  modes  in  the  material  gain  re¬ 
gion  that  spans  from  the  normalized  frequency  of  0.28  to 
0.34,  their  irreducible  representations  in  the  C6v  point 
symmetry  group  are  labeled  in  Fig.  1.  Even  though  this 
paper  focuses  on  the  application  of  group  theory  to  iden¬ 
tify  modes  calculated  from  FDTD,  we  shall  point  out  that 
application  of  the  group  theory  can  also  be  used  during 
the  numerical  simulation  to  restrict  the  resonant  modes 
only  to  those  satisfying  the  symmetry  condition,  thereby 
improving  calculation  convergence. 
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4.  RESONANCE  MODES’  FAR-FIELD 
RADIATION 


The  polarization  of  a  resonant  mode  can  also  be  experi¬ 
mentally  measured  and  numerically  evaluated.  Compar¬ 
ing  the  polarization  measurement  with  the  calculation  re¬ 
sult  serves  as  additional  confirmation  of  the  numerical 
predictions. 

As  shown  in  Figs.  2  and  3,  the  near  field  of  the  lasing 
mode  assumes  a  60-deg  rotational  symmetry.  It  is  there- 
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Fig.  3.  Electrical-field  components  Ex  (top  left)  and  Ey  (top 
right)  at  the  midplane  for  the  19-missing-hole  suspended- 
membrane  photonic-crystal  defect-cavity  mode  considered  (top), 
with  a  two-dimensional  dielectric  contour  overlaid.  Their  projec¬ 
tions  onto  the  irreducible  representations  of  the  C6„  group  are  la¬ 
beled  by  Aj,  A2,  Blt  B2,  Elt  and  E2 ■ 
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Fig.  4.  Magnitude  of  the  Hz  component  at  the  midplane  for  the 
19-missing-hole  suspended-membrane  photonic-crystal  defect- 
cavity  mode  considered  (top),  with  a  two-dimensional  dielectric 
contour  overlaid.  Its  projections  onto  the  irreducible  representa¬ 
tions  of  the  C6„  group  are  labeled  by  A!,  A2,  Blt  B2,  Elt  and  E2. 


fore  expected  that  the  polarization  of  the  radiation  for  the 
lasing  mode  is  not  linearly  polarized.  Its  polarization  of 
the  radiation  varies  little  with  angle;  the  experimental 
characterization  of  the  lasing  mode  is  consistent  with  this 
prediction.  We  have  also  obtained  agreement  between  the 
predicted  spatial  profile  and  the  radiation  pattern  for  this 
mode.14  More-specific  agreement  between  prediction  and 
polarization  measurement  occurs  for  modes  at  the  nor¬ 
malized  frequency  of  0.309  (X0 ~  1618  nm).  Figure  4 
shows  the  magnetic-field  component  Hz  at  the  midplane 
of  the  microcavity.  Its  projections  on  {A1,A2,B1,B2, 
Ei,E2}  of  the  CR,,  point  group  are  shown  in  the  same  lin¬ 
ear  gray  scale.  As  with  the  lasing  mode,  the  resonance  af- 
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ter  FDTD  and  the  Blackman  filter  has  a  negligible  projec¬ 
tion  onto  representations  other  than  E1. 

As  shown  in  Table  1,  E1  of  the  C6v  group  is  a  two- 
dimensional  representation,  which  implies  that  the  eigen¬ 
mode  is  doubly  degenerate.  For  a  two-dimensional  repre¬ 
sentation,  the  final  mode’s  spatial  distribution  depends  on 
the  coefficients  of  its  two-dimensional  linear  combination. 
These  coefficients  depend  on  the  initial  condition.  Figure 
5  shows  another  simulation  of  the  same  Ei  mode  excited 
with  a  different  initial  condition.  Even  though  the  field 
distributions  in  Figs.  4  and  5  are  visually  different,  they 
are  invariant  after  the  same  set  of  symmetry  operations. 


Fig.  5.  Magnitude  of  the  magnetic  field  Hz  at  the  midplane  (top) 
of  the  19-missing-hole  suspended-membrane  photonic-crystal  de¬ 
fect  cavity  for  the  same  mode  as  in  Fig.  5  but  excited  with  differ¬ 
ent  initial  condition.  Its  projections  onto  the  irreducible  repre¬ 
sentations  of  the  C6„  group  are  labeled  by  A1;  A2,  B1,  B2,  Elt  and 
E2' 


Table  2.  Character  Table  for  the  C2v 
Point-Symmetry  Group 

C2„  (2 mm) 

E 

c2 

A” 

1 

1 

1 

1 

A-2 

1 

1 

-1 

-1 

B\ 

1 

-1 

1 

-1 

B’l 

1 

-1 

-1 

1 

The  Ei  mode  is  usually  observed  experimentally  as  two 
closely  spaced  peaks.  We  attribute  the  degeneracy  split¬ 
ting  to  a  perturbation  resulting  from  the  fabrication  pro¬ 
cesses.  The  photonic-crystal  patterns  were  created  by 
electron  beam  lithography  working  in  Cartesian  coordi¬ 
nates.  The  lithography  creates  a  rectangular  perturbation 
that  will  lower  the  symmetry  and  result  in  the  decompo¬ 
sition  of  El  into  the  irreducible  representations  of  the  C2v 
point  group.  The  C2v  point  group  consists  of  four  one¬ 
dimensional  irreducible  representations  A",  A2,  B"2,  and 
B"2,  as  shown  in  Table  2. 

The  fields  after  this  decomposition  are  shown  in  Fig.  6. 
The  E1  mode  of  the  C6v  group  has  zero  projections  on  A" 
and  A',[.  Its  projections  on  B'[  and  B"2  are  identical  regard¬ 
less  of  the  initial  excitation,  which  is  what  is  expected 
from  the  relations  between  the  Ei  of  the  C6v  group  and 
the  B'[  and  the  B2  of  the  C6u  group.  Similarly,  the  double 
degeneracy  of  the  E2  mode  of  the  C6v  group  is  expected  to 
be  split  on  to  the  Aj  and  A2  irreducible  representations  of 
the  C6v  group.  At  the  normalized  frequency  near  0.33, 
there  are  two  doubly  degenerate  modes  labeled  as  E1  and 
E2  in  the  Fig.  1.  Even  though  both  resonances  at  the  fre¬ 
quency  are  displaced  from  the  prediction  owing  to  the  fab¬ 
rication  variation,  they  are  nonetheless  observed  experi¬ 
mentally  as  two  closely  spaced  peaks. 

The  far-held  radiation  in  the  free  space  can  be  consid¬ 
ered  as  an  electromagnetic  problem  with  a  current  or 
charge  distribution  as  a  source.  It  is  often  approached  by 
dyadic  Green’s  function,15 


E(r)  =  icofj,  dVj(r')-G(r,r'),  (1.9) 


where  the  dyadic  Green’s  function  is 


K 


Fig.  6.  Decomposition  of  the  magnetic  field  Hz  as  shown  in  Figs. 
5  and  6,  a  two-dimensional  E1  mode  of  the  C6„  group,  in  one¬ 
dimensional  representation  {.£?'■[, B2}  of  the  C2„  group. 
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Fig.  7.  Calculated  Ex  and  Ey  far-field  components  of  the  modes 
corresponding  to  the  B'[  and  B"2  representations  shown  in  a  polar 
plot.  The  dashed  circle  indicates  a  numerical  aperture  of  0.55, 
identical  to  the  lens  system  employed  for  experimental  light 
collection. 


G(r,r') 


exp(i&0|r  -  r'|) 


(1.10) 


According  to  the  uniqueness  theorem,  a  field  in  a  lossless 
region  is  uniquely  specified  by  the  source  in  it  plus  the 
tangential  components  of  the  electric  field  over  the  bound¬ 
ary,  or  the  tangential  components  of  the  magnetic  field 
over  the  boundary,  or  the  tangential  components  of  the 
electric  field  over  part  of  the  boundary  and  the  tangential 
components  of  the  magnetic  field  over  the  rest  of  the 
boundary.15  Therefore,  the  radiation  from  a  locally  dis¬ 
tributed  source  outside  the  source  can  be  calculated  by 
equivalent  surface  current  or  charge. 

E(r)  =  dVV  X  G(r,r')  •  h  X  2E(r').  (1.11) 

J  S 

By  use  of  the  far-field  approximation, 

|/+ly  VJ  ~(/-ff),  (1.12) 


the  far-field  radiation  of  each  nondegenerate  mode  is  then 
calculated  with  a  vector  Green’s  function  propagation  of 
the  electric  field  Es(r')  at  the  surface  of  the  resonant  cav¬ 
ity: 


ik  o  exp(i£0r) 

E(r)  = - - - r0  X 

4ttt 


-2zX  E.(r') 


Xexp(- ik0  •  r')d2r',  (1.13) 

in  which  k0  is  the  free-space  wave  vector  from  the  center 
of  the  defect. 

Figure  7  shows  the  polar  plots  of  the  Ex  and  Ev  far-field 
components  for  these  two  modes.  The  field  magnitude  at 


the  center  indicates  the  field  radiation  along  the  optical 
axis.  The  field  at  the  same  radiation  angle  is  shown  in  the 
plot  as  concentric  circles.  The  dashed  circle  indicates  the 
numerical  aperture  of  the  lens  being  used  in  the  experi¬ 
mental  system.  Light  emitted  from  the  photonic-crystal 
microcavity  is  vertically  collected  by  an  optical  system 
with  a  numerical  aperture  of  0.55.  For  the  B'[  mode,  the 
Ex  field  goes  through  zero  at  the  optical  axis  and  the 
peaks  of  the  field  locate  outside  of  the  numerical  aperture. 
The  Ey  field,  however,  peaks  around  the  optical  axis.  Be¬ 
cause  the  optical  system  can  collect  only  a  small  angle  of 
the  light  efficiently  from  the  cavity,  the  mode  appears  pre¬ 
dominantly  y  polarized  for  the  lens  system.  In  contrast, 
B"2  is  predominantly  x  polarized. 

This  prediction  of  the  degenerate  E1  modes  splitting 
into  two  modes  that  are  linearly  polarized  along  orthogo¬ 
nal  directions  was  verified  experimentally.  Figure  8  shows 
the  measured  spectrum  and  the  intensity  of  these  cavity 
modes  through  a  polarizer  as  a  function  of  angle.  The  in¬ 
set  shows  the  calculated  far-field  radiation  pattern  of  the 
respective  modes.  It  is  worth  pointing  out  that  the  shapes 
of  the  far-field  intensity  are  not  the  reason  for  determina¬ 
tion  of  the  polarization  of  the  corresponding  mode.  The 
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Fig.  8.  (a)  Measured  optical  spectrum  of  the  cavity,  (b)  Polariza¬ 

tions  of  the  two  resonances  split  in  frequency  during  to  a  rectan¬ 
gular  perturbation  of  a  doubly  degenerate  mode.  The  insets  show 
the  calculated  far-field  intensities  of  the  respective  modes. 
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transmitted  optical  power  peaks  as  the  polarizer  is 
aligned  at  170°  for  the  ffy  mode  and  approximately  60°  for 
the  B2  mode;  they  are  nearly  orthogonal.  Both  modes 
show  a  typical  sine-squared  dependency  on  the  angle  as 
expected  for  linear  polarization.  The  peak-to-valley  ratio 
of  the  optical  intensity  as  the  polarizer  was  rotated  is 
about  10  for  both  cases. 


5.  SUMMARY 

In  summary,  the  resonant  frequencies  and  the  quality  fac¬ 
tors  of  defect  modes  in  a  19-missing-hole  suspended- 
membrane  photonic-crystal  microcavity  are  calculated 
with  3-D  FDTD  and  Pade’s  approximation.  Good  agree¬ 
ment  is  obtained  between  the  calculated  and  the  mea¬ 
sured  resonant  frequencies.  The  resonant  modes  are  clas¬ 
sified  by  the  irreducible  representations  of  the  C6„  point 
group.  The  free-space  far-held  radiation  patterns  are  cal¬ 
culated  with  the  vector  Green’s  function  and  are  com¬ 
pared  favorably  with  the  polarization-resolved  measure¬ 
ments.  Although  this  paper  has  focused  on  a  particular 
19-missing-hole  triangular-lattice  photonic-crystal  cavity, 
the  method  detailed  here  can  be  generally  applied  to  the 
analysis  of  other  multimode  photonic-crystal  microcavi¬ 
ties.  Even  though  the  group  theory  is  used  to  identify 
modes  calculated  from  FDTD  in  this  paper,  it  can  also  be 
applied  in  reverse.  Enforcement  of  the  symmetry  proper¬ 
ties  of  an  irreducible  representation  during  the  numerical 
simulation  restricts  the  resonant  modes  only  to  those  sat¬ 
isfying  the  symmetry  condition,  thus  allowing  a  wider  nu¬ 
merical  filter  to  be  employed  as  the  spacing  between 
modes  is  generally  increased,  and  therefore  accelerating 
the  convergence  to  the  mode  of  interest. 
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